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December 15, 1892. 

Sir JOHN EVANS, K.C.B., D.C.L., LL.D., Treasurer, in the Chair. 

A List of the Presents received was laid on the table, and thanks 
ordered for them. 

The Bight Hon. John Morley, a member of Her Majesty’s Most 
Honourable Privy Council, whose certificate had been suspended as 
required by the Statutes, was balloted for and elected a Fellow of the 
Society. 

The following Papers were read :—• 

I. 44 On an Apparatus for facilitating the Reduction of Tidal 
Observations.” By Gc. H. Darwin, F.R.S., Plumian Pro¬ 
fessor and Fellow of Trinity College, Cambridge*. Received 
November 12, 1892. 

§ 1. Introduction .. 

The tidal oscillation, of the ocean may be represented as the sum of 
a number of simple harmonic waves which go through their periods 
approximately once, twice,, thrice, four times in a mean solar day. 
But these simple harmonic waves may be regarded as being rigorously 
diurnal, semi-diurnal, ter-diurnal, and so forth, if the length of the 
day referred to be adapted to suit the particular wave under con¬ 
sideration. The idea of a series of special scales of time is thus in¬ 
troduced, each time-scale being, appropriate to a special tide. For 
example, the mean interval between successive culminations of the 
moon is 24 h 50™, and this interval may be described as the mean 
lunar day. Now there is a series of tides, bearing the initials M 1? 
M 2 , M 3 , M 4j &c., which go through their periods rigorously once, 
twice, thrice, four times, &c., in a mean lunar day. The solar tides, 
S, proceed according to mean solar time; but, besides mean lunar 
and mean solar times, there are special time scales appropriate to the 
larger (N) and smaller (L) lunar elliptic tides, to the evectional (r), 
to the diurnal (K x ) and semi-diurnal (K a ) luni-solar tides, to the 
lunar diurnal (0), &c. 

The process of reduction consists of the determination of the mean 
height of the water at each of 24 special hours, and subsequent har¬ 
monic analysis. The means are taken over such periods of time that 
the influence of all the tides governed by other special times is 
eliminated. 
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The process by which the special hourly heights have hitherto 
been obtained is the entry of the heights observed at the mean solar 
hours in a schedule so arranged that each entry falls into a column 
appropriate to the nearest special hour. Schedules of this kind were 
prepared by Mr. Roberts for the Indian Government.* The suc¬ 
cessive rearrangements for each sort of special time were made by 
recopying the whole of the observations time after time into a series 
of appropriate schedules. The mere clerical labour of this work is 
enormous, and great care is required to avoid mistakes. 

All this copying might be avoided if the observed heights were 
written on movable pieces. But a year of observation gives 8,760 
hourly heights, and the orderly sorting and re-sorting of nearly 9,000 
pieces of paper or tablets might prove more laborious and more 
treacherous than recopying the figures. 

It occurred to me, however, that the marshalling of movable 
pieces might be reduced to manageable limits if all the 24 observa¬ 
tions pertaining to a single mean solar day were moved together, for 
the movable pieces would be at once reduced to 365, and each piece 
might be of a size convenient to handle. 

The realisation of this plan affords the subject of this paper, and it 
will appear that not only is all desirable accuracy attainable, but that 
the other requisite of subh a scheme is satisfied, namely, that the 
whole computing apparatus shall serve any number of times and for 
any number of places. 

The first idea which naturally occurred was to have narrow sliding 
tablets which should be thrown into their places by a number of 
templates. It is unnecessary to recount all my trials and failures, 
but it will suffice to say that the slides and templates require the 
precision of a mathematical instrument if they are to work satis¬ 
factorily, and that the manufacture would be so expensive as to make 
the price of the instrument prohibitive. 

The idea of making the tablets or strips to slide into their places 
was then abandoned, and the strips are now made with short pins on 
their under sides, so that they can be stuck on to a drawing board iu 
any desired position. The templates, which were also troublesome to 
make, are replaced by large sheets of paper with numbered marks on 

* An edition of these computation forms was reprinted by aid of a grant from 
the Royal Society, and is* sold by the Cambridge Scientific Instrument Company, 
but only about a dozen copies now remain. In the course of the preparation of 
the “ guide sheets ” of the method proposed in this paper, I found that there are 
many small mistakes in these Indian farms, but they are fortunately not of such a 
kind as to produce a sensible vitiation of results. I learn that the mistakes arose 
from a misunderstanding on the part.of a computer employed to draw up the forms. 

The accuracy of my guide sheets was controlled by aid of Mr. Roberts’s forms, 
and it was the occasional discrepancy between my results and the forms which led 
to the detection of the errors referred to. 
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them to show how the strips are to be set.. The guide sheet is laid on 
a drawing hoard, and the pins on the strips pierce the paper and fix 
them in their proper positions. 

The shifting of the strips from one arrangement to the next is cer¬ 
tainly slower than when they slid into their places automatically, 
hut I find that even without practice it only takes about 7 or 8 
minutes to shift 74 of them from any one arrangement to a new 
one. 

The strip belonging to each mean solar day is divided by black 
lines into 24 equal spaces, intended for the entry of the hourly 
heights of water. The strip is 9 in. long by 4 in. wide and the 
divisions (f by 4 ) are of convenient size for the entries. There was 
much difficulty in discovering a good material, but after various 
trials artificial ivory, or xylonite, was found to serve the purpose. 
Xylonite is white, will take writing with Indian ink or pencil, and 
can easily be cleaned with a damp cloth. It is just as easy to write 
with liquid Indian ink as with ordinary ink, which must not be used, 
because it stains the surface. 

The strips have a great tendency to warp, and I have two methods 
of overcoming this. A veneer of xylonite on hard wood serves well, 
or solid xylonite may be stiffened by sheet brass let into a slot on the 
under side. In the first plan the pins are fixed in the wood, and in 
the second the brass is filed to a spike at each end. Whichever plan is 
adopted, the strips are expensive, costing about £7 for a set, and I do 
not at present see any way of making them cheaper. 

The observations are to be treated in groups of two and a half 
lunations or 74 days. A set of strips, therefore, consists of 74, num¬ 
bered from 0 to 73 in small figures on their flat ends. 

If a set be pinned horizontally on a drawing board in vertical 
column, we have a form consisting of rows for each mean solar day 
and columns for each hour. The observed heights of the water are 
then written on the strips. 

When the 24 columns are summed and divided by the number of 
entries we obtain the mean solar hourly mean heights. The har¬ 
monic analysis of these means gives the mean solar tides. But for 
evaluating the other tides the strips must be rearranged, and to this 
point we turn our attention. 

Let us consider a special case, that of mean lunar time. A mean 
lunar hour is about l h 2 m m.s. time ; hence the 12 h of each m.s. day 
must lie within 31 m m.s. time of a mean lunar hour. The following 
sample gives the incidence to the nearest lunar hour of the first few 
days in a year:— 
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The successive 12 h of m.s. time will march retrogressively through 
all the 24 hours of m. lunar time. 

Now, if starting from strip 0, we push strip 1 one division to the 
left, strip 2 two divisions to the left, and so on, the entries on the 
strips will be arranged in columns of approximately lunar time. 

The rule for this arrangement is given by marks on a sheet of 
paper 18 in. broad; these marks consist of parallel numbered steps 
or zigzags showing where the ends of each strip are to be placed so 
as to bring the hourly values into their proper places. 

At the end of a lunation mean solar time has gained a whole day 
over mean lunar time and the 12 h solar again agrees with the 12 h 
lunar. On the guide sheet we see that the zigzag which takes its 
origin at the left end of strip 0 has descended diagonally from right 
to left until it has reached the left margin of the paper, and a new 
zigzag then begins on the right margin. 

When the strips are pinned out following the zigzags on the sheet 
marked M, the entries are arranged in 48 columns, but the number of 
entries in each column is different. The 48 columns are to be 

regarded as appertaining to 0 h , 1\.. 22 h , 23 h , 0 h , l h , . . . ., 22 h , 23 h . 

Thus, the number of entries in the left-hand column of any hour 
added to the number of entries in the right-hand column of the same 
hour is, in each case, 74. The 48 incomplete columns may, in fact, 
be regarded as 24 complete' ones. 

The 24 complete columns are-then summed; the 24 sums would, if 
divided by the total number of strips, give the 24 mean lunar hourly 
heights. The harmonic analysis of these sums gives certain constants, 
which, when divided by the- number of strips, are the required tidal 
constants. It must be remarked, however, that, as the incidence of 
the entries is not exact in lunar time, investigation must be made of 
the corrections arising out of this inexactness. 

The explanation of the guide sheet for lunar time will serve, mutatis 
mutandis , for all the others. 

The zigzags have to be placed so as to bring the columns into exact 
alignment, and printers’ types provide all the accuracy requisite. 
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Accordingly, the computing strips are made to suit a chosen type. 
The standard length for one of the 24 divisions on the strips was 
chosen as that of a “ 2-em English quadrate ”; 24 of these come to 
9 inches, which is the length of a strip. I found the English quadrate 
a little too narrow, and accordingly between each line of quadrates 
there is a “ blind rule,” of 42 to the inch. The depth of the guide sheet 
is that of 74 quadrates and 74 rules, making 15-| in. The computing 
strips are ~ in. broad, and 74 of them occupy 14f in. The excess 
of 15 \ above 14f, or t 7 q in., is necessary to permit the easy arrange¬ 
ment of the strips. 

To guard against the risk of the computer accidentally using the 
wrong sheet, the sheets are printed on coloured paper, the sequence of 
colours being that of the rainbow. The sheets for days 0 to 73 are all 
red; those for days 74 to 74 + 73, or 147, are all yellow; those for 
days 148 to 148 + 73, or 221, are green; those for days 222 to 
222 + 73, or 295, are blue; and those for days 296 to 296 + 73, or 369, 
are violet. 

Thus, when the observations for the first 74 days of the year are 
written on the strips all the sheets will be red; the strips will then 
be cleaned, and the observations for the second 74 days written in, 
when all the guide sheets will be yellow, and so on. 

I must now refer to another considerable abridgment of the process 
of harmonic analysis. It is independent of the method of arrange¬ 
ment just sketched. 

In the Indian computation forms the mean solar hourly heights 
have been found for the whole year, and the observations have been 
rearranged for the evaluation of certain other tides governed by a 
time scale which differs but little- from the mean solar scale. I now 
propose to break the mean solar heights into sets of 30 days, and to 
analyse them, and next to harmonically analyse the 12 sets of har¬ 
monic constituents for annual and semi-annual inequalities. By this 
plan the harmonic constants for 11 different tides are obtained by 
one set of additions. In fact, we now get the annual, semi-annual, 
and solar elliptic tides, which formerly demanded much troublesome 
extra computation. A great saving is secured by this alone, and 
the results are in close agreement with those derived from the old 
method. 

The guide sheets marked S and the computation forms are 
arranged so that the observations are broken up into the proper 
groups of 30 days, and they show the computer how to make the sub¬ 
sequent calculations. 

I have also devised an abridged method of evaluating the tides of 
long period MSf, Mf, Mm. The method is less accurate than that 
followed hitherto, but it appears to give fairly good results, and 
reduces the work to very small dimensions. 
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Before entering on the details of my plan it is proper to mention 
that Dr. Borgen has devised and used a method for attaining the 
same end. He has prepared sheets of tracing paper with diagonal 
lines on them, so arranged that when any sheet is laid on the copy of 
the observations written in daily rows and hourly columns, the 
numbers to be summed are found written between a pair of lines. 
This plan is excellent, but I fear that the difficulty of adding correctly 
in diagonal lines is considerable, and the comparative faintness of 
figures seen through tracing paper may be fatiguing to the eyes. Dr. 
Borgen’s plan is simple and inexpensive, and had I not thought that 
the plan now proposed has considerable advantages I should not have 
brought it forward. 

In the investigations which follow the notation of the Report of 
1883 to the British Association on-harmonic analysis is-used without 
further explanation. 

§ 2. Evaluation of A 0 , Sa, Ssa, S x , S 2 , S 4 , S 6 , T, R, K 2 , K i? P. 

The 24 mean solar hourly heights of water are entered in a 
schedule of 24 columns, with one row for each day, extending to n 
days ; the 24 columns are summed, and the sums divided by n ; the 
24 means are harmonically analysed; it is required to find from the 
results the values of the harmonic constituents. 

The speed of any one of the tides differs from a multiple of 15° per 
hour by a small angle; thus, any one of the tides is expressible in the 
form H cos [(I5°g—/3)£-—£], where q is 0, 1, 2, 3, &e., and j3 is small. 

When t lies between 0 h , and 24 h this formula expresses the oscilla¬ 
tion of level due to this tide on the day 0 of the series of days. 

If multiples of 24 h from 1 to n—l be added to t, the expression gives 
the height at the same hour, t , of mean solar time on each of the suc¬ 
cession of days. 

Then if ^ denotes the mean height of water, as due to this tide 
alone, at the hour t , we have 

1 n 

$ = Voos [(15 °q—p) t — § - +(l5°g-p) 24 (n l)] 
n x 

= -Z cos [(15 .'q-flt-Z-Vtfl (»-l)]. (1). 

n sm p 

When t is put successively equal to 0 h , l h , . . . . , 23 h we get the 24 
values of 1) which are to be submitted to harmonic analysis. 

The mean value of f), say A 0 (not to be confused with A 0 as written 
at the head of this section, where is denoted the mean sea level above 
datum) is found by taking the mean of the 24 values of I). 

By the formula for the summation of a series of cosines it is easy to 
prove that 
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4s ^4S MCt+{iHfl+w - (2) - 

We will now find the jp th harmonic constituents A p , By the 
ordinary rules 


i}=*!0 5 > 


(3). 


Now 


2 V 1 15 °pt COS [(1 Bq-fi) *-^-12 /? («- — l)l 
COS 


= [( 15 °(2+i ? )-/ 3 ) t-Z-12p(n- 1)] 


+ COS 


-sin CO« # (2“*)-/») i-r-12/3 (»-!)] ; 


and T V of the sum of the 24 values corresponding* to £ = 0, 1 ,...., 23 is 

i 8ml2[l5 0 ( ? +p)-^]co8. ris , o f , V 3) _ t B( y, 

2 4 sin£ [15° ( 3+i >)-/ij sin La C 12) ‘ } ? M )J 

+ the same with sign of p charged 

This expression admits of simplification, because 12 x 15° = 180° ; 
making this simplification, and introducing the result into (3), we 
obtain 


B t 


} = — 
J 247 


„. „ J+»“ s k+( 1! >”-4)^+¥”( S +p): 

H sin 12 nB< _!--- 

1 | sin [-V-° (2+p)-i^] 


-dn [?+(“»-*) fl+V (q-p)l 

sin [Y“ (/l-p)-iPi 


.... (4). 


In the particular case where p == q, we have 


A 
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Q 



-Usin 12nB< 

24n 


Ainti+( 12n ~i)P + 15 °‘l] 
Mirpi 1 V=p) 


-j- cos 
+ sin 


[f+(l2«-i) /3] 


( 5 ). 
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If the number of days n be large, Ap, will be small unless the 
denominator of one of the two terms in (4) be very small. This last 
case can only occur wlien p = q and when ft is small. Hence, in the 
analysis of a term of the form under consideration, we may neglect 
all the harmonics except the £ th one. Accordingly (2) and (5) are the 
only formula© required. 

A case, however, which there will be occasion to use hereafter is 
when n = 30, q = 2, when (4) becomes 


A 2 

B 2 


} 


f cos 


, H sin 360 ft 


i' 


(£ + 3591 / 3 ) 
sin |/3 


+ ^Gr+3° 0 + 359i/3)1 


sin (30°—^ ft) 


| ( 6 )- 


For the present we have to apply (5) in the two cases q = 1, 
ft = 0°*0410686 and q = 2, ft = 0 C *0821372 ; now the ratios of cosec f ft 
to cosec (15 °q — \ft) in these two cases are 722 to 1 and 697 to 1. In 
both cases the first term of (5) is negligible compared with the 
second. 


How write 


4 = 


24 n sin \ft 
sin 12 nft 


(7), 


and (5) becomes, with sufficient exactness, 


B, 


} = 


H cos 


[£ + (l2»-i>/3] 


( 8 ). 


If this be compared with (2),. we see that when q = o this formula 
also comprises (2). 

In the applications to be made ft is very small, so that is approxi¬ 
mately a function of the form 6 cosec 0, This function increases very 
rapidly when 6 passes 90°, but for considerable values less than 90° 
it only slightly exceeds unity; for' example, when 0 — 60°, = 1*2, 

but when 6 = 180°, = infin. 

It follows, therefore, that if the number n of days in the series is 
such that 12 nft is less than say '60°, the magnitudes of A q , B q are but 
little diminished by division by ; but if 12^/2 is nearly 180°, A q , B q 
become vanishingly small. 

If the typical tide here considered be the principal lunar tide M 2 , 
and if the number of days be as nearly as possible an exact multiple 
of a semi-lunation, 12 nft is nearly 180°, and the corresponding A 2 , B 2 
become very small. JSTo number of whole days can be an exacl mul¬ 
tiple of a semi-lunation, so that A 2 , B 2 corresponding to M 2 cannot be 
made to vanish completely. For the present they may be treated as 
negligible, and we return to this point in the next section. 

The above investigation shows that in the expression for the whole 
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oscillation of sea level upon which the proposed analysis is performed 
all those tides may be omitted from which ft is not very small, and 
also all those whose frequencies are such that the period under con¬ 
sideration 12nfi is nearly 180°. 

Since the period under consideration will be a lunation, it follows 
that, as far as is now material, the general expression for sea level 
may be written as follows, t denoting mean solar hour angle equal to 
15 °t: — 


m. w., annual ... A 0 +II sa cos (h—ic sa ) 

semi-annual. -f 0^« cos (2 h — K ssa ) 

Solar tides, S a , S 2 . + 0^ cos (t— xf) -f H s cos (2t—**) 

S 3 , S4 . “f" COS (4 t - K2s) “f -035 COS (6 1- K^s) 

Solar elliptic, T . +Ht cos (2t—fr+p,—«$) 

B. +-Hi* COS (2t-f-&— jp t + 7T — K r ) 

Luni-solar, K 3 . -j-f'IT" cos (2t + 2h~2v" — k") 

Ki +f R cos (t + h — v r — g - 7r—«/) 

Solar diurnal, P . -f TL V cos (t — h-\- \ *■—"*). (9). 


This includes all the tides whose initials are written at the head of 
this section. 

It is now necessary to break up the year into 12 equidistant luna¬ 
tions of 30 days. This can be done by the omission of 5 days in 
ordinary years, and of 6 days in leap years. 

If the days of the year are numbered 0 to 364 (365 in leap year), 
the twelve months are as follows : — 

0, 0 d to 29 d ; 1, 30 a to 59 d ; omit 60 d ; 2, 61 d to90 d ; 3, 91 d to I20 d ; 
omit 121 d ; 4, 122 a to 151 d ; 5, 152 d to 181 d ; omit 182 d ; 6, 183 d to 
212 d ; 7, 213 d to 242; 8, 243 d to 272 d ; omit 273 d ; 9, 274 d to 303 d ; 
10, 304 d to 333 d ; omit 334 d ; 11, 335 d to 364 d ; in leap year omit 
365 d . 

The increments of sun’s mean longitude from 0 d 0 h of month 0 up 
to 0 h of the day numbered 0 of each group of days or month are as 
follows 

0, 0°; 1, 30°—0 o, 431; 2, 60°*124; 3, 90 o -0 o *306; 4, 120°*249; 
5, 150°—0°*182 ; 6, 180°*373; 7, 210°-0°*057; 8, 240°-0°-488; 
9, 270°*068; 10, 300 o -0°*364; 11, 330°T91. 

Thus if ho be the sun’s mean longitude at 0 d 0 h of month Q, the 
sun’s mean longitude at 0 d 0 h of month t is & 0 -f 30 °t, with sufficient 
approximation. 
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Now let Y with, appropriate suffix denote the initial “equilibrium 
argument ” at 0 d 0 h of month 0, so that 

Vsa = ho, Vssa = 2h 0 , Y* — — h 0 +p „ Y r = ^, + tt, Y" = 2 h 0 ~2v\ 

Y' — h 0 — v—\ 7T, Yj, = — ^ + |tt ; 

then the general expression (9) for the tide in the month r becomes 

cos (?/£ + Y 5a + 30°t— ) + H ssa COS ( 2yt -f V ssa -f 60°t— K ssa ) 

cos (l5°^— -0$ cos (30 0 t—‘K s ) + H 2 s cos (60 0 t—/c 2s ) 

+ H 3s cos (90° t — k 3s ) 

+ HtCOs[(S0°~-'r})t + y t — 30 O r~^Kt] 

-f - -S"/ cos [ (30° -f- f) t d~ Y f T* 30°t—~ 

+ f "H" COS [(30°+ 2 7j) t + Y" + 60°t— /'] 

+ i'H' cos [15° -f ?/) £ + Y' + 30°t — *'] 

q-i/pCOS [(15°— 7j) t-\-Yp —30°r— /q,J ..... (10). 

Each of these terms falls into the type cos [(l5°g—/*)£—-£], and p 
is in every case either ±r),—2rj , or 0. 

Now, when harmonic analysis of the mean of 30 days is carried 
out, coefficients are introduced. 

Write therefore 

$ 24 X 30 sin \ rj ^ 24 X 30 sin rj 

® 1 ~~ si n~360 k] ’ ^ ~~ sin 720 ?/ ' 

With the known value of tj, 

log Jfj = 0-00483, log = 0*01945. 

In applying the method investigated above, it will be observed 
that a term of any frequency 15 °q-~P only contributes to the 
harmonic constituent of order q. 

Then applying our general rule (8) term by term, and observing 
that 359^ = 14° '76, and 719?/ = 29°*53, the result may be written as 
follows:— 
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2^2^ 1 TT COS . R-t COS / TT I OA ° > A O. W/j\ 

§ 2 M j = sin *' +4, sin ^~ V ' + 3 ° 7 + 14 76 ) 

+ J y X Or-Vr- 30 °t 14°'76) 

ISJ-'-SN.< u >- 


g.,«} = 


TT COS 
; oin K & » 


With, the meaning of the term month in the present context, the 
snn has a mean motion of 30° per month, and each of the first five 
g^s and §Ts is a function with ,a constant part and with annual and 
semi-annual inequalities.’ 

When t has successively the 12 values 0, 1, . . . 11, we have 12 

equidistant values of the ^’s.and U’s. These may be harmonically 
analysed for annual and semi-annual inequalities. 

Suppose that the several coefficients to be determined by harmonic 
analysis are defined by the following equations - 

21 0 (t) = -d 0 + Ai cos 30°t -f Bi sin 30°t + -d 2 cos 60°t + B 2 sin 60°t ; 

DO 

d, 

F, 

fi 



C„1 

c 0 

§1« 

Co 

Cl 

& (T) 

’ = K„ 

> + E x 

§2 (T) j 

eo 
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-» 


^4} TV “ 


V COS 30°T + n 


J> sin 30°t + 


E* 

e 2 


>cos 60°T-f- 


J>sin 60°' 
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Then on comparing (12) with (11) we see that:— 
.do — o ) 

Ai\ HsaCOS, y -^ 2 ^ 

14 76}l B 2 j 


^ siu(^- V “- 20 °' 53); 
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Co 
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f • K k 
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El 1 -H-t COS / TT I T A°'ha\ H-T COS / TT 

> = -r- • (ict—vt-r 14 76) . U r — Vr— 14 *76) ; 

Ei J — < g l sm w ' 4- jf 1 sm v ' ’ 


E 2 

e 2 

e 2 

f 2 


+ f"S’" sin 


}=± f 


(/c"—V"—29°'53) .(13). 


From these equations we get 


i(0,-do ■ 
i (ci +Di) _ 

) f'-ET cos , 

1 ” ix Bin^ 

i(0, + d,) ' 

\ H p cos , 

f “IsinW 

i(Ei-f.) 1 
Ke.+FOJ 

Hr cos , 

i(Bi+f.) 1 
!(ex~F x )j 

t = J 1 sin^ 

i(E 2 -f 2 ) 1 

i"H” cos . 

i (e 2 +F 2 ) J 

4* sin ^ 


(14). 


The tidal constants of the several tides enumerated at the head of 
this section are determinable from these equations. 

Our rule is accordingly to analyse in twelve groups of 30 days, 
and then to analyse the resulting harmonic constituents for annual 
and semi-annual inequalities, combining* the final results according to 
the formula just found. 

The edition of “ guide sheets 99 and computation forms which I have 
drawn up are so arranged as to facilitate the whole process and to 
render it quite straightforward. By this single set of additions it is 
thus possible to evaluate eleven tides and mean water. 


§ 3. Clearance of T, B from perturbation by M a . 

The method of the last section was designed to render all the tides 
insensible excepting those enumerated. But M 2 is so much larger 
than any other tide that there is a small residual disturbance which 
ought to he corrected. 

I have made computations, which I do not give, but which show 
that the disturbance of all the harmonic constituents except gji 2 , 
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is insensible. It is required then to determine the correction to be 
applied to ^ 2 , *§ 2 , and thence those for E, F, e, f. 

Suppose that, when time t is counted from 0 d 0 h of month r, the 
M a tide is expressed by M cos [(30°— 0)t— <f], where f3 = 1°*0158958. 

When means taken over 30 days are harmonically analysed the 
formula (6) gives the contributions to ^ 2 , |§ 2 . As it is now required 
to obliterate these contributions, the signs must be changed, and the 
corrections are 

82LW1 „ . JS(r+869^)-^(f+80- + 86^)\ 


= -risin 360 M ' 


sinf /3 


sin (30° —1/3) 


For reasons which will appear below I now write 
£ = r» (T) -0°-5258. 

Then introducing the value of /3 into (15), I find 


W r) \ - __x 
^2 (T) J _ ,2< 


If sin 5° 43' 4 35 


J (U (t) +4°'689) 
sin0' r 30 i ’i8"-6 


+ lin (^ (r> + ^°-689) 
sin 29° 2^52 


Let denote the value of f m < T > at 0 d U h of month 0, and let 
g^h) = g* OT + ph), and let J^ 2 denote a certain factor whose logarithm is 
0*00849, and let M = f H w . 

In the harmonic analysis for the M 2 tide, considered below in § 6, we 
shall have 


, f H m c. 

A 2 = — cos f«, 


i?, = ^sinf„ 

<S2 


Accordingly 


= * [^sin e(T) +^cS e<T> ]- 


These values of Jf C ? S gW (T> must now be introduced into (16), but 
sm 

the algebraic process need not be given in detail. If we write 


i(s+P)l _ i 4 . 

i(R+Q) / 


sin5° 43'35 (4° 4l'-32) 

sinO° ’ 


J(E-Q)/ sin 29° 29' 3l"’4 
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it follows that 


P — 0-01584, Q == 0-00114, R = 0-00147, S = 0-01611. 


Then, when the substitution of the values of M b t m ( T ) is carried 
’ cos s 


out, we find 
r>9Lh) 


}= 


cos eh) 


■P A +Q ft 


+ sin eh) 


r+Q a +p t > 
~~s ~ 42 +Pv 2 


By the definition of eh) it appears that —eh) is the increment 
of twice the mean moon’s hour angle during the time from 0 d 0 h of 
month 0 up to 0 d 0 h of month t, that is to say eh) — —2 (7 —d)t for the 
time specified. The following table gives the values of eh) and of its 
cosine and sine for each month :— 


No. of days 
Month from epoch 0 


(0- 

to epoch r. 

0h). 

cos 0h). 

sin 0h). 

0 

0 

0° 0' 

1-000 

o-ooo 

1 

30 

11 27 

0*980 

0-199 

2 

.61 

47 16 

0-678 

0-735 

3 

91 

58 43 

0-519 

0-855 

4 

122 

tt-85 28 

-0-079 

0-997 

5 

152 

tt-74 1 

-0-275 

0-961 

6 

183 

tt-38 11 

-0*786 

0-618 

7 

213 

7r —26 44 

-0-893 

0-450 

8 

243 

tt-15 18 

-0-965 

0-264 

9 

274 

tt + 20 32 

-0-937 

-0-351 

10 

304 

w + 31 59 

-0-848 

-0-530 

11 

335 

7 r-j-67 49 

-0-378 

-0-926 


If cos eh), sin eh) are regarded as quantities having annual and semi¬ 
annual inequalities, we may write 

cos eh) = ao + oil cos 30° t + / 3i sin 30°t + a 2 cos 60°t -f- ft sin 60°t + .. .. 
sin eh) — ry 0 + 7 X cos 30°t + sin 30°t + 72 cos 60°t + d 2 sin 60°t +. „.. 

On analysing the numerical values of cos eh), sin eh) by the ordinary 
processes, I find, 


But 


§2 


OC 0 = 

: -0*165, 

70 = 

+ 0-273, 

Oil = 

+ 0-626, 

7l = 

-0*500, 

01 = 

: +0*756, 

*1 = 

+ 0-642, 

= 

+ 0-159, 

72 — 

-0-046, 


: +0-199, 

$2 — 

+ 0-166. 

the 

harmonic 

constituents of 

%2 whei 


when analysed for 
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annual and semi-annual inequality were denoted by E 0 , E 1? F 1? E 2 , F 2 , 
and the constituents of §l 2 were denoted by e 0 , e l5 f x , e 2 , f 2 . Hence tbe 
ten corrections to tbe E’s and F’s are (with an easily intelligible 
alternative notation) 

i } 2 — ( — Pato, i, 2 4 " Q70,1,2) A 2 + (Q«o, 1,2 + P70,1,2) R%, 

&F lf 2 — ( — P/ 3 l, 2 + Q^l, 2) -^-2 + (Q/h, 2 "b P^l, 2) R-Zj 

^ e 0, 1, 2 — ( — R*0, 1, 2 — 1,2) -^2 + ( “ Sao, 1, 2 4 “ ^70, 1, 2) p2j 

^fl, 2 = (—R/h, 2~ S&1, 2) A 2 + ( — 2 + P^l, 2) R ?,- 

On substituting tbe numerical values of cc, (3, 7 , < 5 , P, Q, R, S, 
I find 


Coeffit. of A%. 

Coeffit. of J? 2 . 

£E„ — -(-0*0029 

+ 0*0041 

Seo = -0*0042 

+ 0*0031 

^(Ei+fO = -0*0109 

-0*0091 

«i(E!-f0 = +0*0006 

+ 0*0020 

H^ + FO = -0*0020 

+ 0*0000 

alCei-FO = +0*0091 

-0*0108 

^■2 (E 3 +f 2 ) — —0*0028 

-0*0018 

ai(E 2 -f 2 ) = +0*0002 

+ 0*0012 

^■ 2 ( 62 +F 2 ) — —0*0012 

+0*0001 

«|(e 3 -F 3 ) = +0*0017 

-0*0027 


Most of these corrections are negligible, but tbe four which affect 
the solar elliptic tides T, R must be included, because those tides are 
so small that a small error affects them sensibly. Hence we may 
take, with sufficient accuracy, 

£-|-(ei—Fi) “ +0*009 A 2 —0*011 P 2 , ^2 ( e i+Fi) — —0*002 +. 2 , 

^(Ex + fO = —0*01lA 2 —0*009 J? a , a-J(Ex—f x ) = +0*0006 A 2 + 0*002 R 2 

•. (16*), 

where A 2 , R 2 are the components of tbe M 2 derived from tbe re¬ 
duction of that tide by tbe process of § 6 . 

Provision for these corrections is made in tbe computation forms. 

§ 4. _Evaluation of A 0 , Sa, Si, S 2 , S 4 , S 6 , K 2 , Ki, P, when a complete 
year of observation is not available. 

It is now proposed to consider tbe case where tbe period of ob¬ 
servation is as much as six complete months and less than a complete 
year. 


2 b 2 
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T])e method of the last section apparently depends on the com¬ 
pleteness of the year, yet, with certain modifications, it may be 
rendered available for shorter periods. 

We suppose that so much of the year as is available is broken into 
sets of 30 days by the rules of the last section, and that the means 
are harmonically analysed. The results of such harmonic analysis 
for month (t) are given in (11) of § 3, but for the purpose in band 
they now admit of some simplification. It is clear that it is not 
worth while to evaluate the very small solar elliptic tides T and R 
from a short period of observation. If, then, we denote by Pb) the 
ratio of the cube of the sun’s parallax to its mean parallax at the 
middle of the month (t), the first three terms of the third of (11) 

may be included in the expression P^H S k s . The last term of this 

equation really does involve the solar parallax to some extent, and 
we may, with sufficient approximation, write the third pair of equa¬ 
tions 


cos 

|p a ( T )-r-P( T > J s sin 




sin' 


Let us now consider the value of Pb). The longitude of the solar 
perigee is 281° or —79°, and the ratio of the sun’s parallax to its 
mean parallax is approximately 1 + e t cos (Ji + 79°), and the cube 
of that ratio is 1 + 3 e x cos (h + 79°) or 1 + 0*0504 cos (h + 79°) * 
How the sun’s longitude at the middle of month (t), is 
h 0 + 15° + 30°t; hence 


PO) = l-f 0*0504 cos (7^o + 30°t4-94°) 


and 


JL_ 

Pb> 


1 — 0*0504 cos (K 4- 30°t 4- 94°). 


Thus it is easy to compute the values of 1/Pb) for the successive 
months, when we know h 0 the sun’s mean longitude at 0 d 0 h of the 
month 0. 

The semi-annual tide, being usually small, may be neglected in 
these incomplete observations, and the equations (11) now become 

$o (T > = A 0 cos (ic sa -— Vja— 30°t— 14°’76), 


^'1^1 _ tt COS | f 33- COS / t •fji o i A°.Ha\ 

§!«/ ~ ^sinl^+P sin 0 ~ V _3 ° T ~ 14 76 ) 

+ 5 ' q C ° S ( K? -Y^ + 30°T + 14 o -76), 


= *2 
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rr COS 

— • A -*2 Si 

sm 



rr COS 
- Hbs gin ^ 3 Sj 


1 

jpF) 


1—0*0504 cos (ho + 30°t -j- 94°) 

. (17)- 


When the series of successive values of the g/s and ^’s are har¬ 
monically analysed (by processes which we shall consider shortly) the 
several coefficients resulting from such analysis will be defined by 

$o ( t) = A + A cos SO°r+B 1 sin 30°t, 



'Col 
[ C„ j 

C.1 
c». J 

j> cos 1 sin 30°t, 


l_Bol 

■ B.1 

’ ' l 'e 2 J 

cos 60°T-f^ 2 | sin 60°r, 

§ 3 (T)^_p(T)j 

1 e 0 J 

Mean = J. 4 , 


Mean ^/ T ) = P 4> 


Mean g^h) = A 6 , Mean §J 6 h) = B 6 .... (18). 

Then the subsequent procedure as given in (13) and (14) holds 
good, the only difference being that we do not obtain the semi-annual 
and solar elliptic tides. 

We shall now consider the harmonic analysis of an imperfect series 
of values. 

It must be premised that each monthly value of gjt 2 h), §J 2 h) is to be 
divided by its corresponding Ph) before the analysis is made. 

Suppose that Oh) denotes a function which is subject to semi-annual 
inequality, and that 

Oh) = A 0 + A 3 cos 60°t + B 2 sin 60° t 


Then it is clear that 

0<« = Ao+A* 

0 fl) — Aq + ^A 2 + ^\/3.13 2 , 
0< 2 > = A 0 -iA 2 +iysB 2 , 
&c. &c. 


I now define J) 0 > Pi, P 2 thus :— 

D 0 — 0 (0) + 0 {1) 4 * 0 (2) + ...., 

P 1= =O l0) -fiO (1) -iOh)...., 


P 2 = 0 .0 (0 ) +i-/3 0< 1J +£ v/3 0 (2 >.... 
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If there be n equations and if they be treated by the method of 
least squares, we get 

Do — ^A 0 +A 2 (l +J— 2 * * • • ) d-BgCO-f 2 ^3 + 2 \/S ....), 

Di = A 0 (1 + 2 —i' • • • ) +A 2 (1 + 4 + 4 . • • • ) +B 0 (o+lv^ 3 — 4 \/3. • • • )? 
A = A 0 (o+|v^+^v^....)+A 2 ....) 

+ B 0 (o + f+f....).. 

These are the three equations from which A 0 , A 2 , B 2 are to be 
found. 

A schedule is given below for the formation of D 0j A> D 2 > and a 
table of the solutions of these equations according to the number of 
months available. 

Next, suppose CM denotes a function which is subject to annual 
inequality, and that 

CM = A 0 +Ai cos 30 °t+B 1 sin 30°r. 

Then O (0) = Ao+Aj, 

0 (1} = A 0 +-| \/3 Ai +- B 1? 

(j{2) — Ao+^Aj +^\Z3B 15 

&C.y &C. 

In this case the method of least squares gives 

n 0 = c< o >+g {i) +o {2) .... . 

= nA^A x (1 + 2 ♦.. ) +Bj (o+^+^\/3 

A = 0 ( 0 ) +iv/3 0 ^+| 0 < 3) ..., 

= A 0 (i+^-\/3 + i|. • •. )+Ai (l + f + 4 .. • • ) + Bx (0 + ^v/3 + 5 \/ 3 . . .. )y 

ft = o.O« 0 »+iOt 1 >+Jv / 3O (1 > .. 

= A 0 ( 0 +-| + i \/3 .. .. ) + Aj (0 + |-\/3 + 4 a/ 3. . .. ) + Bj (0 + ^ + f.... )* 

Tables are given below for the formation of D 0 , JD U D 2 , and of the 
solutions of the equations according to the number of months avail¬ 
able. 



'Reduction of incomplete series of results. 

Form for finding Do, D l5 D 2 where the monthly values are subject to semi-annual inequality. 
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In both these forms those monthly values which are unknown are to be treated as being zero. For example, if 9 months 

of observation are available, d\ e, f' will be zero. 
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Rule for finding semi-annual inequality from an incomplete series. 


Number of 




months 





available. 

Coeift. of Dq. 

Coefft. of D x . 

Coeift. of _D 2 - 

6 

A„ - 

+ 0-167 




a 3 = 


+0-333 



b 3 = 



+ 0-333 

7 

A 0 = 

+ 0*148 

-0-037 



A, = 

-0-037 

. +0-259 



B 3 = 



+ 0-333 

8 

A 0 = 

+0-136 

-0*045 

-0*026 


a 2 = 

-0-045 

+0-253 

-0-0L9 


B 2 = 

-0-026 

-0-019 

+ 0-275 

9 

A 0 — 

+ 0*123 

-0-027 

-0-047 


A 2 = 

-0-027 

+ 0-228 

+ 0-011 


b 3 = 

-0-047 

+ 0-011 

+ 0-241 

10 

A 0 = 

+ 0-107 


-0-041 


A 2 — 


+ 0-182 



B 3 = 

-0-041 


+ 0-238 


Rule for finding annual inequality from cm incomplete series. 
Number of 


months 


available. 

Coefft of Z> 0 . 

Coefft. of D x . 

Coefft. of Z> 2 ' 

6 

A„ = +0-977 

-0-326 

-1-215 


A, = -0-326 

+ 0*442 

+ 0-405 


B, = -1-215 

+ 0-405 

+ 1-845 

7 

A* = +0-424 


-0-528 


A, — +0*250 




B x = -0-528 


+ 0-990 

8 

A 0 = +0-226 

+ 0-062 

-0-233 


A 1= +0-062 

+ 0-230 

-0-093 


B x = -0-233 

-0-093 

+ 0-552 

9 

A 0 = +0-146 

+ 0-057 

-0-098 


A x = +0*057 

+0-230 

-0-083 


B 3 = -0-098 

-0*083 

+ 0-326 

10 

A 0 = +0-110 

+ 0-036 

-0-036 


A 3 = +0-036 

+ 0-218 

-0-048 


B t = -0-036 

-0-048 

+ 0-218 
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We tlms get the following rule for the evaluation of A 0 , Si, S 2 , 
S 4 , S 6 , K 2 , K 1? P from 6, 7, 8, 9, or 10 months of observation:— 

Proceed as though the year were complete and find the s and 3p\s* 
for as many months as are available. Reduce the ^2 by multiplica¬ 
tion by 1/PO) or 1—0*0504 cos (7fc 0 -i-30 o T-f 94°). 

Analyse for annual inequality , and ^JPW, ^d/Ph) 

for semi-annual inequality according to the rules for reduction of incom¬ 
plete series just given. 

Complete the reduction as in § 3. 

These rules for reduction do not include the case of 11 months, 
nor the case where any month in the series is incomplete (e.g.fii a 
fortnight’s observation were wanting in one of the months), because 
these cases may be treated thus :—the gjt’s and U’s return to the same 
value at the end of a year, and therefore the case of eleven months 
is the same as that of a missing month at any other part of the year. 
In both these cases we may interpolate the missing gj/s and ^’s and 
treat the year as complete. 

If three or more weeks of observation were missing they might fall 
so as to spoil two months, and in this case we should have an in¬ 
complete series. It is then to be recommended that the equations of 
least squares be formed and the equations solved. So many similar 
cases may arise that it does not seem worth while to solve the equations 
until the case arises. 


§ 5. Evaluation of A 0 , S 2 , S 4 , K 2 , K 1? P from a short period of 
observation. 

If the available tidal observations only extend over a few months, 
it is useless to attempt the independent evaluation of those tides 
which we have hitherto found by means of annual and semi-annual in¬ 
equalities in the monthly harmonic constants. We will suppose that 
30 days of observations are available. Then when we neglect the 
annual tide, and the solar (meteorological) tide Si, we have from (11) 
or (17), which give the analysis of 30 days, 


^0 — A 0 , 



- c ? s ( k '_V'-14°-70)+-J °? S («j-V» + 14°-76), 
sm v y cfi sm V P V J 

= PH- Ks + (k" V" *9°'53), 


IPs gin K„, P = 1+ 0'0504 cos (ho +15°). 


It is now necessary to assume that the P tide has the same 
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amount of retardation as the K 1} and that the ratio of their ampli¬ 
tudes is the -same as in the equilibrium theory. We also make the 
like assumption with respect to the K 2 and S 2 tides. 

Accordingly we put 

Up — <3 H , Kp = K J II ' ~ “k K — &S‘ 

Now since 


Y' = JlQ — \7T—v\ Yp = —/lo + l^r, Y" = 2/?o — 2i/', 

we have 

/^-Y^-fift 6 = /c f —Y' — 14°*76+(2/i 0 — i/'+29 0 -53)+7r, 
K "_ Y"~29°*53 = ^ —(2^o—2^'+29 0 *53). 

Therefore 


I 1 } = -A °° S («'—V'—14° - 76) 
Si / A sm v ' 


i-ETcos 


• (k —y' —14°-76 -4-2/io-// +29°-53), 


(^2 _ T>TT COS I 1^1 ^ ^ CO® /" ^1, [ o 0 /\0.r<v\ 

SJ - P ^ sin + sin ( K « 2^o + 2 1 / -29 53). 

T , , , sin (2h 0 — i/ + 29°-53) 

Let us put tan 0 = — -i—^---i— . 

3f'— cos (2h 0 —/ + 29°-53) 


Then 


tan 0- = _ f> Sln C 2 ^- 2 *" +29°-53) 

' -V--PJf 2 + f" COS (2 h t) — 2 y” +29°'53) 


1 1 _ H’ [3f-C0S (2feo-r' + 29°-53] C OS , , Ar , 


(19). 


3 COS ( 


( K '_V' —14°-76—0), 


-E 

§li - g 

_jr V-Pg + t" COS (2fe„-2 +29°'53) C OS , . s , om 

§J VA'cob* -sin^-^) •• O0). 


If therefore 


l 1 _ T) oos <> 

J - Rl sin ?1 ’ 




we have 


k — £*i-j-Y / + 14 o, 76-f-0 — Kp, 
3 J^Ra cos 0 


H' = 0 -ir 


3f —COS (2fe 0 —v'+29°-53) ’ 


TT — XJT' 

, — 3 • 



367 


1892.] facilitating the Reduction of Tidal Observations , 

KS = £ 2 +^ — lC "> 


T r _ _ -V-JaRaCOS yjr _ 

Us ~~ —P^ + f^COS ( 2 /^ 0 —2i/ / -f-29°*53) ’ 


ff'=*2Z,..(21). 


If there he several months available it is recommended that each 
30 days be treated quite independently, so that from each group of 
days we shall get JET, k' and H s , k s . Then the mean value of 
H' cos k is to be taken as the final value of that function, and H 1 sin k 
is to be treated similarly ; finally ZT, k are to be found. The several 
values of H Si *s may be treated in the same way. Of course we 
assume throughout that k p = P, H p = k" = %, H" = t 8 t-27, as¬ 
sumptions which are usually nearly correct. 

The mean value of must be taken as giving A 0 , but at places 
with a considerable annual tide it is impossible to obtain a good 
value of mean water mark from a short series of observations. 


§ 6. On the evaluation of the several tides by grouping of mean 
solar days. 

Let n( 7~x) denote the speed in degrees per m.s. hour of any one 
tide, n being equal to 1, 2, 3, 4, 5, or 6. Then 15°/(7—%) may be 
called one “special hour.” Since 157(7 — 7 ) is one m.s. hour, the 
ratio of the m.s. to the special hour is ( 7 —x)/(7— 7 ). 

Let one m.s. hour be equal to 1 —/3 special hour, then 

ft = special hours. 

Let it be required to express the 12 h of any m.s. day of a series of 
days by reference to special time. It is clear that 12 h m.s. time will 
be specified by one of the 24 special hours, with something less than 
half a special hour added or subtracted. 

Having fixed the 12 h of m.s. time of a particular m.s. day in 
the special time scale, let us treat that m.s. day as a whole, and 
consider the incidence of the other 23 m.s. hours in special time. It 
is clear that in m.s. time we work backwards and forwards from 12 h 
by subtracting or adding unity, and that in special time we subtract 
or add 1— ft. 

If 12 h m.s. time be where a. lies between special time, 

the following is a schedule of equivalence :— 
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Mean solar 
time. 

0 11 = 
l h = 
2 h = 

Special time." 
(a^-X2 h ) + (*+12 ft) 
(» h —ll h ) + (*+ll ft) 
(® h — 10 h ) + (« + X0 ft) 

ii" = 

(x h —l h ) + 0 + /3) 

12 h = 

* h -fa. 

13 h = 

O h + l h ) + («—/3) 

22 h = 

(® h + 10 ,1 ) + (a:—10/3) 

II 

& 

CO 

CM 

(» Il + ll h ) + (> —11/3) 


In the column of special time it is supposed that 24 l1 is added or 
subtracted, so that the result is less than 24 h . For example, if x 
is 10, the hour column of special time will run 22 h , 23 h , 0 h , . . . ., 9 11 , 
10 h , ll h , . . . ., 20 h , 21 h . 

If the series of days be long x will have all integral values between 
0 and 23 with equal frequency, and since a has all values between 
-}-■§• and with equal frequency, the excess of the solar hour 

above the nearest exact special hour (which may be called the 
error) will have all its possible values with equal frequency. If the 
mean solar hours be arranged in a schedule of columns headed 
0 h , l h , . . . 23 h of special time, each column will be subject to errors 

which follow the same law of frequency. 
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Let abscissae (fig. 1 ) measured from 0 along A'OA represent 
magnitude of a. 

Since a lies between + -J, the limit of the figure is given by 
OA = OA' = i 

If magnitude of error ( i.e . m.s. — special hour), measured in special 
time, be represented by ordinates, a line BOB' at 45° to AO A' re¬ 
presents all the errors which can arise in the incidence of the m.s. 
12 h in the schedule of special time. 

If a line bV be drawn parallel to and above BB' by a distance /?, 
we have a representation of all the errors of incidence of the m.s. II 11 . 
If a series of equidistant parallel lines be drawn above and below 
BB' until there are 12 above and 11 below, then the errors of all the 
m.s. hours are represented, the top one showing the errors of the 
m.s. 0 h and the bottom one the errors of the m.s. 23 h . 

Any special hour corresponds with equal frequency with each solar 
hour, and hence each.mode of error occurs with equal frequency. 

It is now necessary to consider in how many ways an error of 
given magnitude can occur. If in the figure AM represents an error 
of given magnitude, then wherever MN cuts a diagonal line, it shows 
that an error may arise in one way. 

It is thus clear that there are no + errors greater than ^ +12 / 3 , 
and no — errors greater than ^ + 11 /3, and 

Errors of magnitude. 

•J +12 /3 to ^ +11 may arise in 1 way. 


1+11/8 to £ + 10/3 

>> 

2 ways. 

\ -f-10 /3 to -2 -h 9 f3 


3 ways. 

i-10/3 to £-11/3 

JJ 

23 ways. 

£-11/3 to -(1-12/3) 


24 ways. 

;i-12/3) to -(1-11/3) 

>> 

23 ways. 

.2+9/3) to —(£ + 10/3) 

)> 

2 ways. 

;i + 10/3) to -(1 + 11/3) 

?> 

1 way. 


The frequency of error is represented graphically in fig. 2. The 
slope of the two staircases is drawn at 45°, but any other slope would 
have done equally well. 

A frequency curve of this form is not very convenient, and, as 
there are many steps in the ascending and descending slopes, I sub¬ 
stitute the frequency curve shown in fig. 3. This is clearly equivalent 
to the former one. In fig. 3 all the times shown in fig. 2 are con¬ 
verted to angle at 15° to the hour; e accordingly denotes 15°/3. 
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Now let cos n($—x) be the observed value of a function whose 
true value is cos n 0 , and suppose that x, the error of 0, has a fre¬ 
quency / (x) • then the mean value of the function deduced from many 
observations will be 


~+CO p+OO 

/ (x) cosn (0— x) dx-r- /(a?) 
J —oo j —00 


dx. 


In our case /(a?) is the ordinate of the frequency curve whose 
abscissa is x. 

Let OQ = h, QB = a, QB' = b, OA = a + h, OA' = b-\-h ; then 

p + oo 

/ (x) dx = (a + b + fb) h. 

J — CO 


' + 00 

/ (A) cos n (0—x) dx 


i a p«+& 

h cos n (0—x) dx+ (a + h—x) cos n (0—x) dx 

o * a 

rh p5+& 

+ h cos n (5 + 7*.-— x) cosn(0-j-x) dx, 

Jo j 1 } 

4 

= —j cos n [d—§ (a —5)] sin nh sin J n (a + b + h ). 
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The algebraical steps involved in the evaluation of these four 
integrals and subsequent simplification are omitted. 

Hence tbe result is 


sin f nh sin f n (a b + hi) 
f nh ^n^a + b + h) 


cos n (0—\ ( a — b)). 


By reference to the figure it is clear that 
« + & + & = 15°, h ~ 24 e, a =7f°— life, b = ° — 12fe, a—b — c. 

Write, then, 

4 __ 12 ^e - 1 #n 

7 sin 12 we sin-—?* ’ 


and we obtain as the mean value of cos nO , when found in this way, 


4- cosn(0— |e). 

< 31 » 

It is obvious that if we had begun with sin n0, the argument in the 
result and the factor would have been tbe same. Accordingly, a func- 

R/ 

tion R' cos (n0— £') would yield tbe result — cos e) — £']. 

$ n 

If 24 equidistant resulfs of tbis sort are submitted to harmonic 
analysis to find A n , B w , we shall get— 

B/ 

A n == — cos ff' + f ne) — B cos suppose, 

B f 

B n = sin (£■'+f we) = B sin suppose. 

R' , . 

Accordingly R = ^-> ? = ? +i we - 


But it is required to find B/, f', so tbat 

B' = JvjR, ~ ?—f We. 

Thus when tbe 24 observed hourly tide heights on any m.s. day are 
regrouped so tbat tbe observed height at 12 h m.s. time is reputed to 
appertain to an exact special hour, and each of the previous and sub¬ 
sequent hourly values of that m.s. day are reputed to belong to previous 
and subsequent exact special hours ; and when a long series of m.s. 
days are treated similarly, and when the mean height^ of water at 
each of the 24 special hours are harmonically analysed, we shall obtain 
the required result by augmenting B by a factor and by subtract¬ 
ing f ne from £\ 
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The values of and of \nz will be different for each kind of tide, 
and the following table gives their numerical values. 


Table of and \ne. 


Initial of tide. 

n , 

log Jf „. 



1 

0-00212 

0°-26 

m 2 

2 

0-00849 

0°-53 

m 3 

3 

0-01915 

o 

o 

*<r 

to 

m 4 

4 

0*03416 

1°*05 

m 6 

6 

0*07767 

l°-57 

1ST 

2 

0-01361 

GO 

o 

o 

L 

2 

0-00570 

Cq 

0 

O 

V 

2 

0-01278 

0°-78 

0 

1 

0-00535 

0°-57 

J 

1 

0-00225 

—-0 o, 28 

Q 

1 

0-01149 

o 

Oi 

6 

o 

fX 

2 

0-02016 

l°-09 

2SM 

2 

0-00805 

— 0°*49 

MS 

4 

0-02342 

0°-52 

X 

2 

0-00595 

0°*28 

21ST 

2 

0-02136 

1°-13 

00 

1 

0-00481 

—0°*53 

MK 

3 

0-01438 

0°-50 

2MK 

3 

002632 

1°*09 

MN 

4 

0-04328 

l°-35 


As it does not appear worth while to evaluate the tides written 
below the line, no use will be made of the last six results given in 
this table. 

§ 7. On the periods over ivhich the means are to he taiten in evaluating 
the tidal constants. 

We have considered in previous sections the treatment of the group 
of tides which are associated with solar time, when the period of ob¬ 
servation is less than a year, and we have now to consider the other 
tides. 

It is important that the means be taken over such a number of 
days that the perturbation arising from other tides shall be mini¬ 
mised. 

The perturbation between semi-diurnal and diurnal tides is always 
negligible. It is therefore only necessary to consider the action of the 
tides M 2 , S 2 in the case of semi-diurnal tides, and that of Ki and 0 for 
diurnal tides. 
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It is easy to see that the influence of a disturbing tide is evanescent 
when the means are taken over a period such that the excess of the 
argument of the disturbed over that of the disturbing tide has in¬ 
creased through a multiple of 360°. As, however, we are working 
with integral numbers of days, and as the speeds of tides are incom¬ 
mensurable, this condition cannot be exactly satisfied. 

From this consideration it appears that to minimise the perturbation 
of So, 2SM, fx by M 2 (and vice versa ) we must stop at an exact multiple 
of a semi-lunation. To minimise the effect of M 2 on 1ST and L, and of 
K x on J and Q, we must stop at an exact multiple of a lunar anom¬ 
alistic period. To minimise the effect of M 2 on v , we must stop at a 
multiple of the period 2 7 r/((?-\-w — 2 f). To minimise the effect of 
on O, we must stop at an exact multiple of a semi-lunar period. 

For the quater-diurnal tide, MS, it is immaterial where we stop, 
and so it may as well be taken at a multiple of a semi-lunation. 

The following table (p. 374) gives the rules derived from these 
considerations. 


§ 8. On the tides of long period. 

The annual (Sa) and semi-annual (Ssa) tides are evaluated in the 
course of the work by which other important tides are found. These 
are the only two tides of long period which have a practical import¬ 
ance in respect to tidal prediction, but the luni-solar fortnightly 
(MSf), the lunar fortnightly (Mf), and the lunar monthly (Mm) tides 
have a theoretical interest. 

It will therefore be well to show how they may be found. The 
process is short, and, although it is less accurate than the laborious 
plan followed in the Indian reductions, it appears to give fairly good 
results. 

For the sake of simplicity, let us consider the tide MSf. Its 
period is about 14 days, and therefore a day does not differ very 
largely from a twelfth part of the period. Accordingly, if about 
two days in a fortnight are rejected by proper rules, the mean heights 
of water on the remaining days may be taken as representatives of 
twelve equidistant values of water height. 

I therefore go through the whole year and reject, according to 
proper rules, the daily sums of the 24 hourly heights corresponding 
to certain 69 of the days out of 369. The remaining 300 values are 
written consecutively into a schedule of 12 columns and 25 rows, of 
which each corresponds to a half lunation. The 12 columns are 
summed, and the sums are harmonically analysed for the first pair of 
harmonic components. These components have to be divided by 24 
times 25, or by 600, because the daily mean water height is of 
the daily sum, and there are 25 semi-lunations. 

In the same way the semi-lunar period is about 13y days, and if 

VOL. lii. 2 c 
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Number of the last day to be included in the evaluation of the several 
tides for observations extending over any period up to a year . 

For M 2 , ju, 2SM, MS. For O. For N, L, J, Q. For r. 

Stop with one Stop with one Stop with one Stop with one 

of the following of the following of the following of the following 

days days days days 

(semi-lunations), (semi-lunar periods), (anom. periods). (periods 

27rl(ff + m—2rj )). 


14 

13 

29 

26 

43 

40 

58 

54 

73 

67 

74+14 

74+ 7 

+ 28 

+ 21 

+43 

+34 

+ 58 

+48 

+ 73 

+62 

148 + 13 

148+ 1 

+28 

+ 15 

+43 

+ 29 

+ 58 

+42 

+ 72 

+ 56 

— 

+ 70 

222 + 13 

— 

+ 28 

222+ 9 

+43 

+ 23 

+ 58 

+37 

+ 72 

+ 50 


+ 64 

296 + 13 

• - 

+28 

296+ 4 

+ 43 

4-1* 

+ 57 

-r32 

+ 72 

+45 

—__—- 

+ 58 


+ 72 


27 

31 

54 

63 

74+ 8 

74 + 20 

+ 35 

+ 52 

+ 63 

148 + 10 

148 + 16 

+42 

+44 

--- 

+ 71 

222+ 0 


+ 31 

222 + 25 

+ 63 

+ 53 

296 + 21 

296+ 6 

+ 53 

+34 


+ 61 



we erase by proper rules 45 daily sums out of 369, we are left with 
324, which may be written consecutively in a schedule of 12 columns 
and 27 rows, of which each corresponds to a semi-lunar period. The 
summing and analysis is the same as in the last case, but the final 
division is by 24 times 27, or by 648. 
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In this way we evaluate the luni-solar fortnightly and lunar fort¬ 
nightly inequalities in the height of the water. 

The period of the moon is between 27 and 28 days, and if we erase 
appropriately about one day in eight we are left with sets of 24 values 
which may be taken as 24 equidistant values of the daily sums. 
Accordingly we erase 46 daily sums out of 358, and write the 312 
which remain consecutively into a schedule of 24 columns and 13 
rows, of which each corresponds to a lunar anomalistic period. 

The 24 columns are summed and the sums analysed for the first com¬ 
ponents. Finally, the components are to be divided by 24 times 13, or 
by 312. In this way the lunar monthly tide is evaluated. 

But the result obtained in this way is, as far as concerns the tide 
MSf, to some, and it may be to a large, extent fictitious. It repre¬ 
sents, in fact, a residuum of the principal lunar tide M 2 . That this is 
the case will now be proved. 

Suppose that t 0 is an integral number of days since epoch, being 
the time of noon on a certain day ; then the principal lunar tide M 2 on 
that day may be written H m cos [2 ( 7 —<0(^<>+ T ) — fj, where t is less 
than 24 hours. Then the daily sum for that day will be 

-“COS [2 (7 — ff) 4+23 (7—(7) — £»]. 

sin (^ 7 — cr) 

Now since t 0 is an integral number of days 2 (7—0-) t 0 only differs 
from —2 (<r— 7 ) t 0 by an exact multiple of 360° ; hence the argument 
of the cosine may be written 2 {a— 7 ) t 0 — 23 (7— o-) -f £ m . 

But the true luni-solar fortnightly tide, which we may denote 
fflcos [2 (a-— 7 X 4 +t) — £], varies so slowly in the course of a day 
that the daily sum is sensibly equal to 

24 f H cos [2 (o-—7) t 0 +23 (<7—7)—g]. 

It thus appears that the residual effect of M 2 is of exactly the same 
form as that of MSf. It becomes, therefore, necessary to clear the 
harmonic components, determined as described above, from the effects 
of M 2 . 

In order to determine the values of these clearances, I found the 
values of cos 2 (<r— 7 )^ and sin 2 (<r— y)t for every noon in a year of 
369 days. I then erased the values selected for the treatment of 
MSf and analysed the remaining values. In this way it was easy to 
find the effect of the known M 2 tide. 

Suppose that Ai, B x are the first harmonic components determined 
by the treatment of a series of daily sums, and that SAi, ^B x are the 
corrections to be applied to them to eliminate the effects of M 2 , then 
I find that if A m , B m are the two components of M 2 as determined by 
the previous method (§ 6 ) of analysis, 


2 c 2 
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SAi = +0*0304 A m —0*0171 B m , 

= —0*0171 A . m — 0*0304 B„„ 

0 = A x + SA 1} D = B] + £1+. 

C—0*047 D = 0*992 iH cos g*, 

D + 0*047 C — 0*992 fU sin g\ 

Whence f being known from Baird’s manual (being a function of 
the longitude of moon’s node), H and g' are determinable. We have 
also 

tf = f+2(s 0 —If— &o + +11°*7. 

In the set of computation forms which I have prepared for use on 
the present plan, it is shown what days are to be erased for each of 
the three analyses, and how they are to be entered in schedules, 
summed, and analysed. 

§ 9. On abridgment in the computations. 

It seemed probable that one decimal of a foot would suffice to 
express the hourly tide heights. In order to test this, I have taken 
several individual days of observation at Port Blair, and have found, 
by harmonic analysis, the time and amplitude of the diurnal and 
semi-diurnal H.W., first, when the hourly heights are expressed to 
two decimal places of a foot, and secondly, when they are only 
entered to the nearest tenth of a foot. I find that the times of H.W. 
agree within less than a minute of time, and that the amplitudes 
agree within a fraction of an inch. If this much be true of indi¬ 
vidual days, the difference of results arising from two or one place of 
decimals will clearly entirely disappear when a series of days is con¬ 
sidered. Hence, by taking as unit the tenth of a foot, or the inch, or 
even two inches at places with large tides, we may always express all, 
or nearly all, the heights on which we are to operate by two signifi¬ 
cant figures. The adoption of this rule not only saves the writing of 
a large number of figures, but also enormously diminishes the labour 
of the additions w r hich have to be made. 

It also seemed probable that substantial accuracy might be attained 
from the harmonic analysis of only 12 hourly values instead of 24. 
In order to test this I took the tidal reductions for Port Blair, 
Andaman Islands (kindly lent me by the Survey of India), and have 
compared the results which would have been derived from 12 values 
with those actually obtained from 24 values by the computers of the 
Indian Survey. The following tables give the results 
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Semi-diurnal tides . 


Initial. 

' 

Results from 12 
two-hourly values. 

Results from 24 
hourly values. 

Error, (12) —(24). 

H 

A 2 — 

ft. 

+ 0*6883 

ft. 

+ 0-6890 

ft. 

-0*0007 

LB 3 = 

+ 0-6775 

+ 0-6768 

+ 0-0008 

M- 

r A 2 = 

-1-7032 

-1-7005 

-0-0027 

b 2 — 

+1 -0883 

+1 -0872 

+ 0*0011 

rr\ 

rA 2 = 

-0-1437 

-0-1407 

-0 -0030 


b 2 — 

-0-2527 

-0-2515 

-0-0012 

T, 

a 2 = 

+ 0-0357 

+ 0 -0347 

+ 0 -0010 

Ju *< 

LB, = 

+ 0-0612 

+ 0 0610 

+ 0-0002 


rA 2 = 

+ 0*3422 

+ 0-3486 

-0*0064 

LB, - 

+ 0 -2192 

+ 0-2124 

+ 0-0068 

H 

fA 2 = 

-0-1217 

-0*1165 

-0-0052 

LB, = 

+ 0 -0867 

+ 0-0887 

-0-0020 

H 

r a 2 = 

+ 0-0857 

+ 0-0849 

-0*0008 

[B 2 = 

+ 0-0383 

-0-0388 

-0-0005 

2SM J 

rA 2 = 

+ 0-0055 

+ 0-0037 

+ 0*0018 

IB* = 

-0-0198 

-0*0200 

+ 0-0002 


Diurnal tides. 


Initial. 

Result from 12 
two-hourly values. 

Result from 24 
hourly values. 

% 

Error, (12) — (24). 


ft. 

+ 0-0175 
+ 0*0223 

ft. 

+ 0*0185 
+ 0-0216 

ft. 

-0-0010 
+ 0-0007 

K 

H li 

+ 0-0120 
—0 *0168 

+ 0*0059 
-0-0173 

+ 0*0061 
+ 0 -0005 

W 

bdj> 

11 1! 

+ 0-3815 
+ 0-1398 

+ 0*3847 
+ 0-1396 

-0-0032 
+ 0-0002 

O 

11 11 

-0-0818 
+ 0*1335 

-0-0729 
+ 0-1386 

-0-0089 

-0-0051 


-0-0167 

-0-1287 

-0*0178 

-0-1280 

+ 0*0011 
-0-0007 


-0*0193 
+ 0-0315 

-0-0167 
+ 0-0347 

-0-0026 

-0-0032 

«{£:: 

+ 0*0140 
-0 -0170 

+ 0-0136 
-0-0194 

+ 0-0004 
+ 0-0024 
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The mean discrepancy in the case of the semi-diurnal tides is 
0*0022 ft., and the greatest is +0*0068; in the case of the diurnal 
tides the mean discrepancy is 0*0026 ft., and the greatest is 0*0089. 

In tidal work results derived from different years of observation 
differ far more than do these two sets of results, and hence the ana¬ 
lysis of 12 two-hourly values for diurnal and semi-diurnal tides gives 
adequate results. 

I find that this abbreviation does not give satisfactory results for 
quater-diurnal tides, and the sixth harmonic is not derivable from 12 
values. Therefore, when these tides are to be evaluated the 24 
hourly values must be used. 

It will still be necessary to write all the 24 hourly heights on each 
computing strip, but when the strips are put into any one of the 
arrangements, except where quater-diurnal tides are required, we 
need only add up the columns 0, 2, 4, ...., 22, and may omit the 
columns 1, 3, ...., 23. 


§ 10. On a trial of the proposed method of reduction. 

As already mentioned, I have the tidal reductions for one year 
(beginning April 19, 1880) for Port Blair, Andaman Islands. I am 
thus able to make a comparison between the results of the old method 
and of the new. The computation was, in large part, done for me by 
Mr. Wright. 

It appeared sufficient to evaluate the tides of the S series and 
those allied with them, the tides of the M series, and the tide Q; also 
the tides of long period MSf, Mf, Mm. 

The S series test the new process of harmonic analysis of monthly 
harmonic components for annual and semi-annual inequalities. I 
chose M because it is the most important tide, and Q because it puts 
the proposed method of grouping to a severe test, and is very small in 
amplitude. 

In the Q time scale the day is 26 h 52 m of mean solar time, from 
which it follows that one of the 24 mean solar hourly observations 
may fall as much as 2 h 0 m away from the exact Q hour to which it is 
reputed to belong. Thus the hourly observations are arranged in 
wide groups round the Q hours, and the hypothesis involved in the 
method is put to a severe strain. 

Lastly, the results for tides of long period test my proposed 
abridgment. 

It will be seen in the table on p. 379 that the two methods give 
results in close agreement. There is, however, a sensible discrepancy 
in the K 2 tide, but in this case I am inclined to accept the new value 
as better than the old one. This tide is governed by sidereal time, 
which differs but little from mean solar time. Hence, in the Indian 
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Tort Blair . 1880-81. 




I. 

Indian 

calculation. 

1ST. 

New 

method. 

I-N 

(height). 

I-N 
(phase). 



ft. 

ft. 

ft. 


■A-o 

= 

4 -792 

4-795 

-0-003 


*{? 

_ 

*299 

•299 

000 


= 

163° 

162° 


+ 1° 

« fH 

_ 

*106 

Til 

- -005 


tbsa ■< 

= 

165° 

164° 


+ 1° 

*{? 

_ 

•099* 

•094 



= 

313°* 

339° 





•020* 

•004 




326°* 

312° 



*{? 


•028 

•026 

+ -002 


= 

49° 

53° 

•• 

-4° 

*{? 

= 

*966 

•973 

- -007 


= 

316° 

315° 

*• 

+ 1° 



•003 

*003 

000 


— 

107° 

105° 

... 

+ 2° 

Mf 

= 

•403 

•401. 

+ *002 


= 

326° 

326° 

•• 

0° 

Mf 

— 

•286 

*268 

+ -018 


— 

314° 

311° 

•• 

+ 3° 

pjs 

— 

•130 

•139 

- -009 



= 

324° 

323° 

•• 

+ 1° 


_ 

•014 

•013 

+ -001 


- 

23° 

34° 

• • 

-11° 


= 

2-042 

2*043 

- *001 


= 

279° 

279° 

• • 

0 

m 3 {h 

= 

*004 

•004 

000 


= 

20° 

54° 

i 

-34° 


_ 

•003 

*006 

- -003 


= 

167° 

264° 

j 

-97° 


_ 

•004 

•005 

- *001 


= 

342° 

315° 

•• 

+ 27° 

(24 values) Q-j^ 

= 

*023 

236° 

*023 

233° 

000 

i- 3° 

(12 values) Q, ^ 

= 

•023 

236° 

*022 

234° 

+ *001 

+ 2° 


* These are derived from 1880-82. 
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Tides of Long Period. 



I. 

Indian 

calculation. 

N. 

New 

method. 

I-N 
(height). 

I-N 
(phase). 


ft. 

ft. 

ft. 


= 

•045 

163° 

•019 

168° 

+ *026 

-5° 

K 

* M 

it it 

•056 

*056 

000 


356° 

356° 

1 

0 

Lk = 

*016 

12° 

•020 

13° 

- *004 

-1° 


method of grouping, considerable errors of incidence of the S hours 
in the K time scale prevail for many days together, and the method 
seems of doubtful propriety. The same is true of the P tide, and here 
also the two methods give somewhat different resuits. 

The accuracy with which the very small Q tide comes out, whether 
from 24 values or only from 12, is surprising, and may perhaps be, to 
some extent, due to accident. It shows, however, that the present 
method may be safely applied, even when the special time scale 
differs considerably from mean solar time. 

The results for the tides of long period are quite as close to the old 
values as could be expected. 

§ 11.-4 comparison of the ivork involved in the new and old methods of 

reduction. 

It has been usual in the Indian reductions to use three digits in 
expressing the height of water, and there have been 15 series, or even 
more. Now 3 x 24 x 365 x 15 is 394000 ; hence the computer has had 
to write that number of figures in reducing a year of observation. 
This does not include the evaluation of the annual and semi-annual 
tides, so that we may say that there have been about 400,000 figures 
to write. 

I propose to express the heights by two digits, and they only have 
to be written once. Thus, in the present plan, the number of figures 
to write is 2x24x365, or 17,500. Thus the writing of 382,000 
figures is saved. 

In the old method the computer had to add together all the digits 
written, say, 394,000 additions of digit to digit. 

I propose to use 24 hourly values in three series, viz., S, M, and 
MS, and 12 two-hourly values in eight others. Therefore, the number 
of additions will be 3 x 2 x 24 x 365 + 8 x 2 x 12 x 365 or 123,000. Thus 
270,000 additions are saved. 
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We may say that formerly there were about 800,000 operations 
(writing and addition), and that in the present method there will be 
about 140,000. This estimate does not include a saving of several 
thousands of operations in obtaining the tides of long period. I am 
therefore within the mark when I claim that the work formerly be¬ 
stowed on one year of observation will now reduce at least five 
years. 

It has been found that the manufacture of my computing strips of 
xylonite is rather expensive, but as it formerly cost in England 
rather more than d£20 to reduce a year of observation, the cost of the 
apparatus will be covered by the saving in the reduction of a single 
year, and it will serve for any length of time. 


§ 12. On the completion of the record for short gaps and long gaps. 

In any long series of tidal observations there are usually some 
breaks in the record in consequence of the stoppage of the clock of 
the tide gauge, or from some other cause. How the process of elimi¬ 
nation by grouping depends essentially on the completeness of the 
record, and it is therefore necessary to fill in blanks by interpo¬ 
lation. 

Such interpolation has not been usual in the operations of the 
Indian Survey, and it might he thought that the complete omission of 
the missing entries is the proper course to take ; but it is easily shown 
that this treatment is exactly equivalent to the assumption that the 
water remained stagnant at mean sea level during the whole time of 
stoppage of the gauge. It is obvious, therefore, that any conjectural 
values are better than none. 

The process by which it is proposed to interpolate is best shown 
by an example. 

At Port Blair (beginning April 19, 1880) the column of 6 11 from 
days 99 to 112 gives the heights shown in the first column of the 
table below. I suppose that the tide gauge broke down on day 103,, 
and only came into action again on day 110.* There was really no 
break down, and the actuality during the supposed hiatus is shown 
in the last column but one. 

How if we look back about a month we find that the water stood 
about the same height at the same hour of the day (viz., 6 h ). Then 
the “ previous record ” (which is complete) beginning at 69 d is en¬ 
tered in the next column. Similarly a u subsequent record ” is found 
about a month later, and is entered in a third column. The mean of 
the previous and subsequent records is then taken as giving the 
values to be interpolated. 

* The days are here numbered from 1, instead of from 0. This has been the 
usage in India hitherto. 
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The last two columns contain a comparison between the interpo¬ 
lation and what in the present case we know to have been actuality. 
There is a mean error of O’20 ft. Thus it is clear that a fair record 
may be interpolated even with so long a break as a week. 

In this example I have only shown the interpolation for one 
column, but of course all the other twenty-three columns would really 
have to be treated similarly. 

I find by trial that the result would be a little improved by a 
graphical method, but that process is slightly more troublesome than 
the numerical one. 


Table of Interpolation. 


Defective 

record. 

Previous 

record. 

Subsequent 

record. 

Mean of 
previous and 
subsequent. 

Actuality. 

Error. 

Day. 

6 h . 

Day. 

6 h . 

Day. 

6-. 





99 

2-54 

69 

2-02 

129 

2’56 

2’29 




100 

3*13 

70 

2-83 

130 

3-07 

2*95 




101 

3*86 

71 

3-70 

131 

3-70 

3-70 




102 

4*40 

72 

4 55 

132 

4-27 

4-41 




103 

,. 

73 

5*10 

133 

4*88 

4-991 

Ji 

4*83 

+ 0*16 

104 

,. 

74 

5-60 

134 

5-27 

5*44 

o 

5*24 

+ 0*20 

105 


75 

5 72 

135 

5*35 

5*54 

-U 

c3 

5*39 

+ 0-15 

106 

,, 

76 

5-67 

136 

5-28 

5’48 

- O 

5*18 

+ 0*30 

107 

• « 

77 

5-59 

137 

4*92 

5*26 

& 

4*90 

+ 0*36 

108 

*. 

78 

5*04 

138 

4 *44 

4*74 

"3 

4*53 

+ 0’21 

109 

,. 

79 

4*62 

139 

3*57 

4 *10 J 

HH 

1 

4*05 

+ 0*05 

110 

3-32 

80 

3-81 

140 

2*95 

3*38 




111 

2*64 

81 

3-26 

141 

2*50 

2*38 




112 

2-17 

82 

2-73 

142 

1*99 

2*36 





It may happen that the hiatus is too long for treatment in this 
way. I do not think it would be safe to treat much more than a fort¬ 
night by interpolation. 

It has been shown in § 4 how the tides associated with S are to be 
treated where the record is deficient, and it remains to consider the 
other tides. 

In § 7 are given the days with which we must stop in the analysis 
of an incomplete year, and this table affords us the means of treating 
a long hiatus in the observation. 

We may in fact omit all the entries between any two of the 
numbers given in the table without seriously affecting the result;. 

Let us suppose, as an example, that the tide gauge broke down on 
day 210 and was only repaired and in operation again on day 226. 
How 210 is 148 + 62, and 225 is 222 + 3. 

Then we see by the table in § 7 that in finding the means for 
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M, 2SM, MS, when the computing strips are written for the third 
time, we must remove strips 59, 60, 61 (which have numbers written 
on them) and may leave the remaining strips of that writing which 
are blank. When the strips are written for the fourth time strips 
0, 1, 2, 3 will be blank, but we mast remove strips 4 to 13 inclusive. 
When all the strips are used in a complete year there are 369, and 
this is the divisor used in obtaining the harmonic constants, but 
when there is this supposed hiatus we do not use 15 strips of the 
third writing and 14 strips of the fourth writing’, so that the divisor 
will be 340. 

Again, when we are evaluating 0 in the third writing, strips 57, 
58, 59, 60, 61 must be removed, and in the fourth writing strips 
4 to 9 inclusive. In a complete year the divisor is 369, but we now 
do not use 17 strips of the third writing and 10 of the fourth writing, 
so that the divisor becomes 342. 

Again, in evaluating N, L, J, Q, in the third writing we remove 
strips 45 to 61 inclusive, and in the fourth writing strips 4 to 25 
inclusive. The divisor is reduced from 358 to 303. 

Lastly in evaluating v, in the third writing strips 43 to 61 inclu¬ 
sive are removed, and in the fourth writing strips 4 to 31 inclusive. 
The divisor is reduced from 350 to 287. 

Any hiatus, be it long or short, may be treated in this way, but it is 
clear that if it be short enough to treat by interpolation, it is best to 
adopt that method. 

INSTRUCTIONS FOR USING THE COMPUTING 
APPARATUS. 

The apparatus for the reduction of tidal observations, together with 
computation forms, can be purchased from the Cambridge Scientific 
Instrument Company at a price (as far as can be now foreseen) of 
about £8. 

In case of any insufficiency in the following instructions recourse 
must be taken to the preceding paper. 

On the degree of accuracy requisite in the hourly heights. 

It will usually be sufficient if the heights be measured to within 
one-tenth of a foot, and the decimal point may, of course, be omitted 
in computation. 

This gives amply sufficient accuracy at a place where the semi- 
range of the principal lunar tide is 2 ft., and where spring range is 
from 6 ft. to 7 ft. 

At some places with small tides a smaller unit might be necessary, 
and at others with very large tides a unit of 2 in., or of a fifth of a 
foot, might suffice. 
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Whatever unit of length be taken it is important, for the saving of 
work, and it is sufficient, that all or nearly all the heights should be 
expressed by two digits. 


Gompletion of record . 

If there is an accidental break in the record, it is very important 
that it should be completed according to the method shown in § 12, 
or by some other equivalent plan. 

The computation forms are drawn up on the supposition that the 
year of observation is complete, but with proper alterations, which 
will now be indicated, they may be used in other cases. 

In §4 it is shown how to treat the tides of the S group when the 
observations have been subject to a long stoppage in the course of the 
year, and also when the observations extend over any period from six 
months to a year. 

In § 5 it is shown how to treat the S group for a short period of 
observation. 

If the stoppage be a long one, the method explained in § 12 must 
be adopted for all the other tides. The same section also shows the 
treatment for observations extending over any period, long or short, 
less than a year. 


Untries and summations. 

The computing strips are intended to take writing in pencil or 
liquid Indian ink , but not in common ink. 

They are to be cleaned with a damp cloth, and a little soda may be 
put in the water if they become greasy. 

Lay the red S sheet on one drawing board and set up the strips 
with their ends abutting against the corresponding numbers. The 
strip numbered 60 is also to be put on the board. 

Write the hourly heights for each day on the strip bearing the 
corresponding number, strip 0 for day 0, strip 1 for day 1, and so on 
up to strip 73 for day 73. The 24 hourly heights are to be written in 
the 24 divisions of each strip, beginning on the left with 0 h and end¬ 
ing on the right with 23 h . 

Remove strip 60. 

Sum the 24 columns formed by the divisional marks on consecutive 
sets of 30 strips. Thus, days 0 to 29 afford 24 sums ; days 30 to 59 
afford the second set of 24 sums ; days 61 to 73 afford 24 sums, which 
are the beginning of a third 30, to be completed when the second set 
of 74 days shall have been written on the strips. 

The numbers 0, 1, 2, ... ., 23, 0,1, . . . ., 23 at the head and foot of 
the guide sheet indicate the hours corresponding to the columns. 

The sums of the columns on the board are to be entered in the 
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corresponding columns of tlie form “ Hourly sums of S series in 
twelve months.” 

Lay the red M guide sheet on the other drawing hoard, and 
transfer the strips from the first hoard to the new arrangement shown 
hy the zigzag lines, strip 60 being now reintroduced. 

There will now be 48 columns (more exactly 47, since one of the 
columns will he found to have nothing in it), numbered at top and 
bottom 0, 1, . . . ., 23, 0, 1, . . . ., 23. Each of the 48 columns is to be 
summed from bottom to top (not as for S in groups of 30), and the 
sums are to be entered in the form “ Sums of series M.” The 24 
sums which come from the left half of the board will be entered in 
the row marked “ red left,” those from the right in the row marked 
“ red right.” 

Lay the red 1ST sheet on the other board, and transfer the strips. 

In accordance with § 9, it will now usually suffice to sum only the 
columns appertaining to the even hours 0, 2, 4, . . . ,, 22; as these 
hours are repeated twice, there will now be 24 columns to sum. 

The sums are then to be entered on the form “ Sums of series H,” 
in the alternate columns. The complete form is provided, so that all 
the 24 hourly values may be used if it be thought desirable, but this 
labour seems unnecessary, at least in a long series of observations. 

Lay the successive red sheets on the vacant board, transfer, sum, and 
enter, until all the red sheets are exhausted. 

In the case of S, M, MS the sums of all the columns are necessary, 
but in the other eight arrangements only the sums of the alternate 
columns, those of the even hours, are usually necessary. For a short 
series of observations it may be best to use all the columns, but in this 
case it will certainly not be worth while to attempt the evaluation of 
v , J, Q, fi, 2SM, which are all small in amount. 

If the tides of long period MSf, Mf, Mm are required, the 24 
numbers written on each strip must be added together, and the sum 
entered in the form “ Long period tides—daily sums.” 

Clean the strips. 

In exactly the same way work through the next 74 days, from 74 d 
to 74 d + 73 d , with yellow guide sheets. Then clean the strips, and 
take another 74 days with green guide sheets, and so on with the 
blue and violet. 

In the last (violet) set attention must be paid to the rules as to the 
places where the analysis is to stop in each arrangement. 

If the year of observation is so incomplete that the hiatus cannot 
be made good by interpolation, or if the series does not run over the 
complete year, the series must stop with one of the days specified in 
the table at the end of § 7, and a note must be made of the number of 
days used in each series. 

The strips marked for omission on the violet sheets, or those 
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selected for omission under the rules of § 7, may be hidden by a sheet 
of paper when the summations are being made. 

The additions of S, M, MS may be verified by proving’ that the 
grand total of all the numbers (inclusive of omitted strips in S) 
written in each of the sets of 74 days is the same in whatever way 
they are arranged. Thus, the sum of the 48 columns should be equal 
to the sum of the daily sums. An incomplete verification in the other 
arrangements, when only half the columns are summed, is found by 
showing that the sum of all the hourly sums of each 74 days is nearly 
equal to half the grand total of all the numbers written in that period 
of 74 days. 

When the guide sheets become worn with many pin pricks they 
may easily be patched with adhesive paper. There seems no reason 
why this patching should not go on almost indefinitely.^ 

* It is possible that it may be desired to evaluate the tides 00 and 2 Y, for which no 
guide sheets are provided. I therefore, give instructions for the preparation of guide 
sheets for these cases. They will be understood by any one who has the set of 
guide sheets before him. With the instructions given below, the computer might 
indeed set up the strips without a guide-sheet. 

I describe the staircase as descending from left to right or from right to left, and 

1 define a short step as being one space down and one space to the left or right, as 
the case may be, and a long step as one space down and two to the left or right, as 
the case may be. When I say, for example, that a short follows 2, I mean that 

2 to 3 is a short step. The first mark on each sheet is specified by its incidence in- 
the row of hours at the top. 

00 ; descending from left to right. 

The sequence is long several times repeated and then short. 

Eed; 0 between O h and l h ; shorts follow 2, 7, 13, 19, 24, 30, 36, 41, 47, 53, 58, 
64, 69. 

Yellow ; 0 between 15 h and 16 h ; shorts follow 1, 7, 12, 18, 24, 29, 35, 41, 46, 52, 
57, 63, 69. 

Green ; 0 between 6 h and 7 h ; shorts follow 0, 6, 12, 17, 23, 29, 34, 40, 45, 51, 57, 
62, 68. 

Blue; 0 between 21 h and 22 h ; shorts follow 0, 5, 11, 17, 22, 28, 33, 39, 45, 50, 
56, 62, 67. 

Violet; 0 between ll h and 12 h ; shorts follow 5, 10, 16, 21, 27, 33, 38, 44, 50, 55* 
61, 67, 72. 

The last strip used for a year is 72. 

2Y : descending from right to left. 

The sequence is long, long, short, long, long, short, and at intervals three longs 
and a short. 

Eed; 0 between 22 h and 23 h ; shorts follow 1, 4, .. .., 16; 20, 23, ...., 35; 39, 
42, ...., 57; 61, 64, 70. 

Yellow ; 0 between 18 h and 19 h ; shorts follow 2 ; 6, 9, ...., 21; 25, 28, ...., 40 j 
44, 47, .. .., 59 ; 63, 66, ...., 72. 

Green ; 0 between 13 h and 14 h ; shorts follow 1,4; 8,11,.. .., 26 ; 30, 33,.. .., 45; 
49, 52, .... , 64; 68, 7l. 

Blue; 0 between 8 h and 9 h ; shorts follow 0, 3, 6, 9 ; 13, 16, .. .., 28; 32* 
35, ...., 50; 54, 57, ....,69; 73. 
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Hourly sums and harmonic analysis. 

Complete tlie summations in the forms for hourly sums, and copy 
into the forms for harmonic analysis. In this copying it will gene¬ 
rally suffice if the last figure in the hourly sums he omitted; for 
example, if the observations are entered to the nearest tenth of a foot 
the hourly sums will be given in the same unit, and it will suffice if 
the hourly sums analysed be written to the nearest foot. 

There are 12 analyses (one for each month of 30 days) for the 
hourly sums in S, and one analysis for each of the other 10 arrange¬ 
ments. All the forms are provided with spaces for 24 hourly sums, but 
in the eight series 1ST, L, v , O, Q, J, ja, 2SM, where only 12 values will 
commonly be used, the entries will only be made on the alternate 
rows of 0 h , 2 h ,. . . 22 h . In these cases the divisor 12, which occurs 

in the penultimate stage of finding the A’s and B’s, must be replaced 
by 6. 

The large divisors (viz., 369 for M, /t, 2SM, MS ; 369 for O; 358 
for N, L, J, Q ; and 350 for v) represent the number of days under 
reduction, and must be altered appropriately (see table, § 7) if there 
be long gaps in the observations, or if the year be incomplete, or if 
the series be a short one. 

If some one of the monthly analyses of S is deficient the missing 
J^’s and §’s are to be made good by interpolation.* 

It is then necessary to analyse the monthly values of the and 
§’s derived from the 12 analyses of S. We thus obtain A 0 , A 1? B x , 
A 2 , B 2 , Co, Co, Cl, Di, Cj, di, Bo, e 0 , E x , B 2 , e x , i u E 2 , E 2 , e 2 , f 2 ,. The rules 
for these analyses when the year is incomplete are given in § 4, and 
the computation forms only apply to the case of the complete year. 

Astronomical data and final reduction . 

Determine from the ‘Nautical Almanac’ and Major Baird’s 
‘Manual of Tidal Observations’* the astronomical data at 0 11 local 
M.T. on day 0, and proceed according to the form to find tne initial 
arguments and factors for reduction. The astronomical data are then 
to be used in the forms for final reduction. 

We have generally B = Using*, A = It cos g*; the forms are 
arranged so that colog A is to be added to log B to find log tan g*, and 
thence g\ If g' lies between —45° and 45° or between 135° and 225°, 
log sec g* is added to log A to find R; if g* lies between 45° and 135° 
or between 225° and 315°, log cosec g’ is added to log B to find R. 

Yiolet; 0 between 4 11 and 5 h ; shorts follow 2, 5,...., 14; 18, 21,.. . •, 33 ; 37, 
40,...., 62 s 56, 59, ...., 71. 

The last strip used for a year is 61. 

* Taylor and Francis, London, 1886, price 7s, Gd, 
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Accordingly the computation form has log sec room being left 
for the syllable “ co ” if necessary; underneath this is written log , 
and the computer will insert A or B as the case may be. 

There is usually required also a numerical factor l/f or Jf, or both ; 
the logarithms of the l/f are found amongst the astronomical data, 
and the logarithms of the constant Jf’s are printed in the forms in 
their proper places. 

The treatment of the 21 harmonic components derived from the 
harmonic analysis of the 5 j^,’s and $J’s is shown in the forms. 

The tides of long period. 

The processes involved in the evaluation of these tides are suffi¬ 
ciently shown in the forms. 


Postscript. December 17, 1892. 

Correction to 'previous paper. 

An error has been detected on p. 333 of my paper “ On the Har¬ 
monic Analysis of Tidal Observations of High and Low Water,” 
‘ Boy. Soc. Proc.,’ vol. 48 (1890). 

In the example of reduction certain tabular values, extracted from 
Baird’s 4 Manual,’ have the wrong signs attributed to them. Since 
in 1887 the longitude of the moon’s node lay between 0° and 180°, 
the signs of r, f, v , 2v" should be positive, instead of negative as 
stated on p. 333. When this error is corrected the reduction leads to 
the following comparison:— 

Computed. Mean of 9 yrs. obs. 
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' P68 

1*63 
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3° 
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= 

0*46 
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= 
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rH 

— 
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1*00 
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= 

317° 

313° 

L i 

” H 

— 

0T1 

0*09 
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= 
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£0 

o 

GO 
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rH 

” 

1*24 

1*40 

l K 

= 

41° 

45° 


'H 

= 

0*69 

0*66 

„ K 

= 

59° 

o 

OO 

P <1 

iH 

= 

0*41 

0*40 
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_ K 

= 

41° 
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If the calculations had been conducted by rigorous methods, the 
two columns would have agreed nearly with one another. 

I may mention that I have copies of a table of ( 7 — 0 -) t up to 
90 days (see p. 304 of the paper here referred to) which I shall be 
glad to give to any one actually engaged in the reduction of H. and 
L.W. observations. 


II. “ On some new Reptiles from the Elgin Sandstone.” By 
E. T. Newton, F.G.S. Communicated by Sir Archibald 
Geikie, F.R.S. Received November 28, 1892. 

(Abstract.) 

During the last few years a number of Reptilian remains have been 
obtained from the Elgin Sandstone at Cuttie’s Hillock, near Elgin, 
which are now in the possession of the Elgin Museum and of the Geologi¬ 
cal Survey. These specimens represent at least eight distinct skeletons, 
seven of which undoubtedly belong to the Dicynodontia, and one is a 
singular horned Reptile, new to science. All the remains yet found 
in this quarry are in the condition of hollow moulds, the bones them¬ 
selves having entirely disappeared. In order, therefore, to render the 
specimens available for study, it was necessary, in the first place, so to 
display and preserve these cavities that casts might be taken which 
would reproduce the form of the original bones. Gutta-percha was 
found to be the most suitable material for taking these impressions; 
and in some instances, especially in the case of the skulls, the casts 
had to be made in several parts and afterwards joined together. 

The first specimen described is named Gordonia Traquairi; it is the 
one noticed by Dr. Traquair in 1885, and referred to the Dicyn¬ 
odontia ; besides the skull, it includes fragmentary portions of other 
parts of the skeleton, and is contained in a block of sandstone which 
has been split open so as to divide the skull almost vertically and 
longitudinally. The two halves have been so developed that casts made 
from them exhibit the left side and upper surface, as well as the main 
parts of the palate and lower jaw. In general appearance this skull 
resembles those of Dicynodon and Oudenodon. The nasal openings are 
double and directed laterally ; the orbits are large and look somewhat 
forwards and upwards. The supra-temporal fossa is large, and 
bounded above by the prominent parieto-squamosal crest, and below 
by the wide supra-temporal bar, which extends downwards posteriorly 
to form the long pedicle for the articulation of the lower jaw. There 
is no lower temporal bar. The maxilla is directed downwards and 
forwards to end in a small tusk. Seen from, above, the skull is 
narrow in the inter-orbital and nasal regions, but wide posteriorly 
across the temporal bars, although the brain-case itself is very 
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